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ABSTRACT 

A  process  model  is  proposed  to  explain  the  various  stages  of  flow  conditions  that  are 
observed  by  measurement  as  the  flow  transitions  from  non-cavitation  to  cavitation 
(turbulent  flow),  supercavitation,  and  finally  separation  in  sharp-edge  90  degree  orifices. 
This  study  includes  orifice  L/D  from  1  to  10,  orifice  diameters  of  nominally  0.048  and 
0.078  inch,  and  upstream  pressures  from  100  to  1500  psi  as  well  as  manifold  cross  flow 
velocity  of  from  6  to  60  ft/sec.  The  results  support  two  different  first  order  models,  one 
for  cavitation  and  the  other  non-cavitation  in  turbulent  flow.  Under  full  cavitation 
conditions  the  discharge  coefficient  is  related  to  the  contraction  coefficient  and  the 
cavitation  number  to  the  V2  power.  In  the  non-cavitation  full  flow  regime  the  head  loss  is 
related  to  the  loss  coefficient  and  the  dynamic  pressure  at  the  orifice  exit.  Both  the  head 
loss  and  contraction  coefficient  were  found  to  be  a  strong  function  of  the  ratio  of 
upstream  velocity-to-orifice  exit  velocity.  The  area  ratio  between  the  manifold  and  the 
orifice  were  also  found  to  have  a  significantly  influence  the  contraction  coefficient. 
Relationships  are  proposed  to  explain  the  processes  leading  to  hydraulic  flip  (separation) 
and  prediction  of  occurrences  that  include  inception  of  cavitation,  supercavitation,  and 
separation. 
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INTRODUCTION 


Rocket  engines  operate  at  both  sea  level  (SL)  and  vacuum  environments.  Depending  on 
mission  requirements  the  engine  may  also  be  required  to  startup  and/or  restart  in  space. 
Conditions  at  both  SL  and  space  result  in  operation  during  startup  within  the  cavitation 
regime  for  injector  designs  incorporating  sharp-edge  orifices.  Orifice  cavitation  can 
cause  damage  to  the  thrust  chamber  and  injector  due  to  unacceptable  jet  characteristics 
(i.e.  misimpingement  of  impinging  jets,  jet  wall  impingement,  and/or  initiation  of 
combustion  instability).  Cavitation  under  these  conditions  can  be  controlled  by  tailoring 
the  startup  sequence  (upstream  pressure/orifice  pressure  drop)  or  design  of  injector 
elements  that  avoid  cavitation  at  all  operating  conditions  (i.e.  rounded  entrances).  It  is 
the  responsibility  of  the  injector  design  engineer  to  ensure  that  the  injector  operates  at 
acceptable  conditions  under  all  operating  conditions.  In  order  to  accomplish  this  a-priori 
rather  than  use  cut-and-try,  the  design  engineer  must  have  design  information  for  (1) 
predicting  manifold/orifice  flow  characteristics  in  both  the  non-cavitation  and  cavitation 
regimes  and  (2)  the  resulting  impact  on  jet  and  flow  characteristics.  To  date  this  type  of 
information  has  been  limited  or  nonexistent. 

Recognizing  this  need  the  US  Air  Force  Research  Propulsion  Laboratory  at  Edwards  Air 
Force  Base,  California  initiated  a  broad  study  of  orifice  flow  in  both  the  cavitation  and 
non-cavitation  regime,  in  a  unique  experimental  facility,  using  a  range  of  impinging 
injector  element  design  variables  consistent  with  rocket  engine  design.  Results  from  this 
facility  were  published  by  Strakey  and  Talley1  in  1999. 

This  paper  describes  the  results  from  90  degree  sharp-edge  orifices  with  manifold  feeder 
cross  flow  from  6  to  60  ft/sec,  L/D  from  1  to  10  and  two  orifice  sizes  (nominally  0.048 
and  0.078  inch).  The  results  of  this  study  provide  an  engineering  process  for  detennining 
the  flow  regime  and  analysis  method  for  defining  the  pertinent  flow  variables. 

COLD  FLOW  INJECTOR  CHARACTERIZATION  FACILITY 
FACILITY  DESIGN 

The  experimental  investigation  was  carried  out  at  the  Air  Force  Research  Laboratory 
cold-flow  injector  characterization  facility,  a  simplified  schematic  of  which  is  shown  in 
Figure  1.  Water,  which  is  used  as  a  simulant  for  liquid  oxygen,  is  stored  and  pressurized 
in  a  35  ft3  tank.  The  injector  inlet  flow  rate  is  controlled  with  a  throttling  valve  and 
measured  with  one  of  several  turbine  flow  meters  arranged  in  parallel  to  cover  a  wide 
range  of  flow  rates.  Downstream  of  the  injector,  another  turbine  flow  meter  measures  the 
outlet  flow  rate  and  a  back-pressure  regulator  maintains  the  injector  fluid  pressure.  Fluid 
pressures  of  up  to  1500  lb/in2  and  manifold  velocities  up  to  130  f/s  may  be  produced  in 
this  manner.  The  injector,  as  illustrated  in  Figure  2,  consists  of  an  interchangeable 
stainless  steel  plate  with  a  sharp-edged  hole  machined  in  the  center  of  the  plate.  This  is 
mounted  against  the  fluid  manifold,  which  contains  a  0.25  in  square  channel  machined 
into  the  manifold.  A  plenum  and  screen  at  the  inlet  and  exit  of  the  manifold  reduces  the 
velocity  before  the  flow  enters  the  manifold  in  order  to  generate  a  reproducible  flow  field 
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at  the  entrance  to  the  orifice.  The  distance  from  the  manifold  inlet  screen  to  the  orifice 
entrance  is  2.5  in  or  9  manifold  channel  widths.  The  entire  arrangement  is  secured  inside 
an  optically  accessible  pressure  chamber  which  is  rated  to  2000  psi. 


Pump 


Figure  1  -  Schematic  of  AFRL  Cold  Flow  Test  Facility 

The  orifices  are  either  pilot  drilled  and  reamed,  or  made  by  electrical  discharge 
machining  with  a  diameter  tolerance  of  +/-  13  pm  and  inlet  edge  radius  to  orifice 
diameter  ratio  of  less  than  .003,  ensuring  a  sharp  edged  inlet.  Chamber  pressure,  orifice 
pressure-drop,  and  inlet  and  outlet  flow  rates  are  recorded  by  a  12  bit  analog  to  digital 
conversion  board  and  the  data  is  stored  on  a  personal  computer.  Experiments  are 
typically  conducted  by  setting  the  fluid  pressure  and  flow  rates  to  a  predetermined  value, 
with  the  chamber  pressure  being  gradually  increased  while  the  data  acquisition  system 
records  flow  rates  and  pressures.  This  allows  for  a  large  amount  of  data  to  be  collected  in 
a  relatively  short  period  of  time.  Orifice  pressure  drop  and  chamber  gas  pressure  are 
measured  within  an  accuracy  of  +/-  0.25%.  Manifold  velocities  are  held  constant  to 
within  +/-  1.5%  during  the  experiment.  Because  of  the  difficulty  associated  with  directly 
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Figure  2  -  Injector  Schematic 

measuring  the  orifice  flow  rate  inside  of  the  pressurized  vessel,  the  orifice  flow  rate  is 
measured  by  subtracting  the  manifold  outlet  flow  rate  from  the  manifold  inlet  flow  rate. 
The  experimental  error  associated  with  the  discharge  coefficient  measurement  is  limited 
by  the  accuracy  of  the  flow  meters  which  is  +/-  0.5%.  This  translates  to  an  error  on  the 
discharge  coefficient  data  of  about  +/-  0.5%  at  the  lowest  manifold  flow  rates  and  highest 
orifice  flow  rates,  to  +/-  10%  at  the  lowest  orifice  flow  rates  and  highest  manifold  flows. 
A  typical  error  for  the  intermediate  flow  rates  is  on  the  order  of  +/-  4%. 

The  range  of  operating  conditions  and  orifice  geometries  typical  of  liquid  rocket  injectors 
was  studied  and  are  given  in  Table  1. 


Table  1  -  Test  Variables  and  Range  of  Testing 


Orifice  Diameter  (in) 

Length  to  Diameter  Ratio 
Fluid  Pressure  (psi) 

Back-Pressure  (psi) 

Cross- Velocity  (f/s) 

Orifice  Red 
Manifold  Rew 

Manifold  Dynamic  Head  to  Orifice  Static  Dp  Ratio 


0.0468  -  .0839 
3  -  10 
100  -  1500 
13.05  -  1500 
0-60.7 

l.OxlO4-  3.0xl05 
6.0xl03 

5.0xl0'5-  6.0x10° 


MANIFOLD  CONFIGURATIONS  TESTED 


Tests  were  conducted  to  distinguish  between  flow  that:  (1)  entered  a  manifold  from  both 
ends  and  then  contracted  to  enter  a  90  degree  sharp-edge  orifice  (Configuration  #  1),  (2) 
entered  the  manifold  from  one  end  and  contracted  to  enter  a  90  degree  sharp  edge  orifice 
(Configuration  #  2),  and  (3)  entered  the  manifold  from  one  end  and  part  of  the  flow 
contracted  to  enter  a  90  degree  sharp  edge  orifice  and  the  remainder  exited  through  the 
manifold  exit  (Configuration  #3).  A  sketch  of  the  configurations  is  shown  in  Figure  3. 
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CONFIGURATIONS 


#1  #2  #3 

Dead  Head  Approach  Velocity  Cross  Velocity 


Orifice  Flow  Out 

Figure  3  -  Orifice  Entrance  Test  Configurations 


For  Configuration  #  1  and  #2  the  ratio  between  the  velocity  V i  of  the  manifold  inlet  and 
the  orifice  exit  velocity  V2  is  constant  and  equal  to  the  ratio  between  the  orifice  area  A2 
and  the  manifold  area  Ai.  This  remains  true  for  Configuration  #3,  but  in  that  case  the 
area  Afis  to  be  interpreted  to  be  the  cross  sectional  area  only  of  the  flow  that  enters  the 
orifice,  not  the  entire  area  of  the  manifold.  This  later  case  therefore  introduces  another 
independent  variable  not  present  in  the  first  two  configurations. 

ANALYSIS  APPROACH 


CAVITATION  TURBULENT  FLOW 


To  account  for  all  the  dynamic  variables  in  the  flow,  the  analysis  below  is  carried  out  in 
terms  of  the  following  dimensionless  variables:  the  area  ratio  A1/A2,  the  velocity  ratios 
V1/V2,  ( V 1  -V 0ut)/V 2  for  Configuration  #3,  where  Vout  is  the  manifold  exit  velocity,  the 
discharge  coefficient  Cd,  the  contraction  coefficient  Cc,  and  the  cavitation  number  Kcav. 
The  discharge  coefficient  is  defined  to  be  the  ratio  of  the  actual  flowrate  to  the  ideal 
flowrate.  The  contraction  coefficient  is  defined  to  be  the  ratio  of  the  area  of  the  vena- 
contracta  Ac  to  the  exit  area  A2.  The  cavitation  number  is  defined  as: 


K 


cav 


Pi~Pr 

P1-P2 


(1) 


Pi  is  the  manifold  pressure,  P2  is  the  backpressure  at  the  orifice  exit,  and  Pv  is  the  liquid 
vapor  pressure.  The  derivation  of  relationships  between  these  dimensionless  variables 
depends  on  various  assumptions  unique  to  the  flow  characteristics  of  the  configurations 
tested. 


In-Line  Orifices 

The  classic  textbook  case  concerns  in-line  orifices  where  Vi  and  V2  both  flow  in  the 
same  direction  and  no  bulk  flow  turning  occurs.  Furthermore,  V 1  is  assumed  small  and 
can  be  neglected.  The  discharge  coefficient  is  then  defined  as: 
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(2) 


Q  =  [te/pX^J'2 

Combining  Equations  1  and  2  results  in: 

c  KJ2Vi 
rf"[(2  glp\Px-PvT 


(3) 


For  this  configuration  the  Bernoulli  equation  together  with  Equations  4  and  5  results  in 
Equation  6: 


C'  =  f- i 

V. 


V.  = 


M 


1/2 


1- 


'V2 


-il/2 


Cc2 


(4) 

(5) 

(6) 


It  should  be  noted  that  the  assumption  made  in  the  application  of  the  Bernoulli  Equation 
is  that  at  cavitation,  the  vena-contracta  static  pressure  is: 

P,  *  P,  (7) 

And  that  the  Hl  term  between  the  manifold  and  the  vena  contracta  is  small. 

Neglecting  the  manifold  velocity  implies  that  the  ratio  A2/A1  is  small  and  that  the  second 
expression  in  the  denominator  of  Equation  (6)  is  small  with  respect  to  1.  Then  for 
engineering  purposes  Equation  6  reduces  to: 

c^cXkJ12  (8) 

2  •  3  10 

This  is  the  expression  used  in  Nurick“  and  subsequent  studies  "  for  In-Line  orifices  and 
found  to  represent  the  data. 

Configuration  #1 

For  all  three  configurations  considered  in  the  present  study  V 1  could  not  be  assumed  to  be 
negligible  and  therefore  it  must  be  included  in  the  determination  of  Cd.  For 
Configuration  #1  and  #2  the  resulting  expression  for  the  discharge  coefficient  is: 
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(9) 


[(2  giplP,-P,)+v;]n 

And  combining  Equations  1,  4,  5,  and  the  Bernoulli  Equation  results  in  the  expression 
showed  in  Equation  10: 


Cd  = 


CK 


1/2 


( ,0 

1  + 

A, 

v  1  y 

1) 


1/2 


(10) 


Note  that  inclusion  of  the  manifold  velocity  accounts  for  the  additional  Kcav  tenn  in 
Equation  10.  Further,  it  was  assumed  in  the  derivation  that  the  major  head  loss  is 
downstream  of  the  vena-contracta,  as  was  assumed  for  the  In-Line  orifices.  For  the 
designs  investigated  in  this  study  the  denominator  for  configuration  #1  and  #2  vary  by 
less  than  1  %  and  therefore  Equation  10  could  be  reduced  to  that  of  Equation  8.  Equation 
10  was,  however,  used  to  define  the  Cd  relationship. 

The  data  shown  in  Figure  4  is  typical  of  the  results  from  Configuration  #  1  testing  using 
Equation  10.  Note  that  a  straight  line  may  be  fitted  through  the  data  in  the  manner  of 
Equation  8  with  a  slope  of  0.6168  and  a  correlation  coefficient  R2  of  0.9999.  Since  the 
head  loss  varies  as  V2  it  would  be  expected  that  if  significant  the  curve  would  not  be 
linear.  The  linear  curve,  therefore,  confirms  that  exclusion  of  the  head  loss  tenn  in 
Equation  10  does  not  impact  the  functional  relationship.  Figure  4  also  supports  that  the 
impact  of  Vi  is  negligible  and  identifies  the  head  loss  post  vent-contracta  as  the  major 
head  loss  area  for  this  configuration. 


Figure  4  -Relationship  between  Cd,  Cc,  and  Kcav  for  Configuration  #1 
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Configuration  #2 

In  this  configuration  all  the  flow  turns  90  deg  from  one  direction.  These  turning  losses 
can  be  significant  if  the  flow  stagnates  against  a  wall  before  turning  into  the  orifice  (i.e. 
for  90  degree  bends  up  to  0.75  times  the  kinetic  head11).  In  Configuration  #2,  where  no 
wall  exists,  the  flow  does  not  stagnate  but  smoothly  turns  to  flow  into  the  orifice.  This 
results  in  a  significant  reduction  in  the  turning  losses  (i.e.  similar  to  a  smooth  bend  ). 
Consequently,  it  was  initially  assumed  that  Equation  10  would  still  be  representative  of 
the  flow  characteristics  for  this  configuration. 

Typical  results  from  Configuration  #2  are  presented  in  Figure  5.  Once  again  a  straight 
line  could  be  plotted  through  the  data  in  the  manner  of  Equation  8  with  a  large  correlation 
coefficient  implying,  as  in  configuration  #1,  that  the  effect  of  V)  as  well  as  the  turning 
losses  was  negligible. 


Figure  5  -  Relationship  between  Cd,  Cc,  and  Kcav  for  Configuration  #2 


Configuration  #3 

For  conditions  where  the  flow  entering  the  manifold  exits  by  two  paths,  (i.e.  the  manifold 
exit  and  the  orifice),  and  the  total  pressure  in  the  manifold  is  constant  the  discharge 
coefficient  can  still  be  defined  by  Equation  9.  However,  Equation  5  then  becomes: 


V„  = 


0 \-vqJa^ 


c. 


v^2  J 


(11) 


Utilizing  the  same  equations  the  resulting  discharge  coefficient  is  now  defined  as: 
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Cd  = 


CK. 


fa  -vout) 


2  {  .  \2 
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vA  j 
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/  \  2 
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n  1/2 


Cc2(KIa,-l) 


The  data  were  first  evaluated  using  Equation  12  to  determine  the  validity  of  the 
assumptions.  If  the  results  were  found  not  to  be  valid  then  the  head  loss  would  be 
included  in  the  analysis.  The  data  shown  in  Figure  6  depict  the  results  using  Equation  12 
for  a  cross  flow  velocity  of  40  ft/sec  and  is  typical  of  the  data  for  other  cross  flow 
velocities. 


Figure  6  -  Relationship  between  Cd,  C„  and  Kcav 

It  was  initially  thought  that  the  turning  losses  would  be  significant  and  therefore  impact 
head  loss.  Any  impact  of  Kcav  and  HL  would  be  revealed  by  a  change  in  the  slope  shown 
in  Figure  6,  suggesting  a  variable  Cc.  The  fact  that  the  slope  is  constant  shows  that  the 
major  head  loss  is  again  downstream  of  the  vena-contracta  and  more  importantly  that  the 
flow  is  choked.  For  all  tests,  choked  flow  is  indicated  by  constant  flowrate  through  the 
orifice  in  the  cavitation  regime.  It  should  be  noted  that  each  test  was  run  at  a  specific 
cross  flow  resulting  in  constant  velocity  ratio. 

For  any  test  where  the  velocity  ratio  remained  constant  there  was  a  unique  value  of  Cc 
(i.e.  as  defined  by  Equation  12).  However,  variations  in  velocity  ratio  should  impact  the 
vena-contracta  formation  process  and  therefore  Cc.  For  Configuration  #3,  V)  can  be 
varied  for  a  given  orifice  flow  by  simply  changing  the  amount  of  flow  that  exits  the 
manifold.  This  results  in  a  significant  change  of  velocity  ratio  (between  0.1  and  1). 
Therefore,  determination  of  Cc  for  Configuration  #3  requires  the  use  on  Equation  12 
rearranged  as: 
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Equation  13  was  used  to  define  Cc  for  all  Configuration  #3  tests. 

NON-CAVITATION  TURBULENT  FLOW 
All  Configurations 

Orifice  discharge  coefficient  in  turbulent  non-cavitation  flow  from  sharp  edge  orifices 
has  traditionally  been  correlated  in  terms  of  Cd  vs  Reynolds  number.  These  correlations 
have  been  developed  for  flow  entering  the  orifice  with  no  cross  flow  and  very  low 
entrance  velocity.  For  configurations  where  the  flow  must  bend  to  enter  the  orifice  or 
contract,  researchers13'14  have  typically  chosen  to  correlate  head  loss  (rather  than  Cd)  with 
the  dynamic  pressure  ( p  V2  /2g).  The  correlating  constant  is  tenned  the  loss  coefficient 
(Kl).  Therefore,  the  total  head  loss  for  turning  constant  velocity  flow  is  defined  by: 

Hl=Kl-?-V?  (14) 

2g 

For  conditions  where  the  velocity  is  not  constant  then  Kl  is  also  be  a  function  of  V1/V2 
or: 

r  y  \ 

H,=K,^-^V2  (15) 

\vj2g 

The  head  loss  can  independently  be  determined  using  the  Bernoulli  equation: 

HL=(Pl-P2)  +  ^(V12-V22)  (16) 

2g 

The  application  of  Equations  15  and  16  to  the  specific  configurations  evaluated  define  the 
Kl  functionality.  In  Configurations  #1  and  #2  the  velocity  ratio  is  constant  and  therefore 
this  variable  drops  out.  For  Configuration  #  3,  Kl  is  a  function  of  velocity  ratio. 

FLOW  CHARACTERISTICS  OF  CONFIGURATIONS  TESTED 

The  various  stages  or  regimes  of  flow  encountered  in  this  study  are  shown  in  Figure 
7.  There  are  5  regimes  and  are  defined  as: 

1.  Non-Cavitation 
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2.  Cavitation  Inception 

3.  Cavitation  with  reattachment  point  moving  downstream 

4.  Supercavitation  with  quasi-reattachment  near  or  at  the  orifice  exit 

5.  Separation 

All  of  these  stages  have  been  noted  by  other  researchers3, 7’ and  8.  In  the  test  indicated  in 
blue  the  flow  is  initially  separated  (Cd  ~  0.6)  then  at  about  a  cavitation  number  of  1.8 
(Kcav°'5  =  1 .35)  flips  to  attached  flow  in  the  non-cavitation  regime.  A  repeat  test  indicated 
in  pink  starts  in  the  supercavitation  regime  (quasi-attached)  until  it  quickly  transitions  to 
full  cavitation,  and  continues  in  cavitation  until  reaching  the  incipient  cavitation  point 
where  the  flow  transitions  to  non-cavitation.  Note  that  2  thru  4  above  all  occur  in  the 
cavitation  regime.  These  characteristics  are  discussed  in  detail  in  the  following  sections. 


Figure  7  -  Results  Showing  the  Various  Stages  of  Flow  Characteristics 
NON-CAVITATION  FLOW  REGIME 

In  the  non-cavitation  regime  the  acceleration  forces  are  insufficient  to  reach  vapor 
pressure  in  the  entry  region  of  the  sharp-edge  orifice.  Nevertheless  the  accelerating  fluid 
in  the  core  separates  from  the  sharp  edge  entrance  producing  a  fluid  recirculation  region 
between  the  separation  and  the  reattachment  point.  The  minimum  static  pressure  occurs 
at  the  vena-contracta.  As  the  acceleration  forces  decrease,  the  core  flow  area  increases 
and  if  sufficient  orifice  length  is  available  reattachment  occurs.  After  reattachment  the 
orifice  pressure  then  decreases  at  a  rate  defined  primarily  by  turbulent  and  friction  losses 
such  that  the  exit  pressure  matches  that  of  the  back  pressure.  The  static  pressure  and 
orifice  flow  characteristics  are  illustrated  below  in  Figure  8. 
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Orifice  Flow  Characteristics 


W1,P1,V1  ' 


Wout,  Pout,  Vout 
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Figure  8  -  Illustration  of  Non-Cavitation  Orifice  Flow  Characteristics 

Typical  results  for  the  90  degree  sharp  edge  orifice  with  Configuration  #1  and  #2  are 
presented  in  Figure  9.  Note  that  the  measured  loss  coefficient  for  both  configurations  is 
only  slightly  different.  The  mean  deviation  of  <  0.001  for  each  curve  supports,  that  for 
these  configurations  (i.e.  no  cross  flows),  the  validity  of  assumptions  in  Equation  10. 


Orifice  Exit  Dynamic  Pressure  ~  psi 


Figure  9  -  Head  Loss  Relationship  in  Non-Cavitation  Turbulent  Flow 

There  is  no  existing  correlation  for  the  90  degree  turning  in  a  cross  flow  manifold  but  it 
would  be  expected  to  be  near  that  for  a  smooth  bend  (i.e.  vanes  consistent  with 
Configuration  2)  which  is  0.2.  Based  on  existing  correlations,  the  head  loss  coefficient 
for  a  contraction  of  0.0769  is  0.48.  For  Configuration  #  2  the  loss  coefficient  is  only 
slightly  higher  (i.e.  0.48+0.2=0.68  vs  0.743)  and  does  not  appear  to  be  inconsistent  with 
expectations.  Interestingly,  KL  is  0.75  for  a  90  degree  elbow  at  constant  velocity11. 
Since  the  only  difference  is  the  contraction,  this  suggests  that  the  major  head  loss  in  non- 
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cavitation  turbulent  flow  is  in  turning  the  flow.  It  follows  then  that  the  major 
contribution  to  head  loss  occurs  upstream  of  the  vena-contracta  since  the  turning  is 
complete  by  that  point  in  the  flow. 

Typical  results  for  the  90  degree  sharp  edge  orifice  with  the  Cross  Flow  Configuration  #3 
are  presented  in  Figure  10.  Each  data  set  was  conducted  at  constant  velocity  ratio  even 
though  the  inlet  manifold  velocity  varied.  Note  that  the  measured  loss  coefficient  varies 
with  manifold  velocity  and  that  increasing  the  velocity  causes  a  greater  deviation  from 
the  zero  intercept.  This  is  a  reflection  of  the  impact  of  VjAC  on  KL.  Further,  for  the 
cross  velocity  configuration  the  loss  coefficient,  Kl,  is  constant  for  any  given  velocity 
ratio  (V1/V2),  and  the  value  of  Kl  varies  with  velocity  ratio.  The  R“  fit  of  the  data  shows 
excellent  linear  correlation.  Consequently,  Kl  will  be  nonlinear  with  respect  to  velocity 
ratio. 


Figure  10  -  Head  Loss  Relationship  in  Non-Cavitation  Turbulent  Flow 

Data  for  the  90  degree  sharp  edge  orifice  where  the  velocity  ratio  varied  from  -0.08  to 
2.5  and  L/D  from  1  to  10  are  shown  in  Figure  11  thru  13.  For  the  Cross  Velocity 
configuration,  90  degree  sharp  edge  orifice,  the  relationship  is: 
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L/D  =  10 
Angle  =  90  deg  I 


KL  =  0.9744(V1/V2)  +  0.3561  (V1/V2)  +  0.6578 
R2  =  0.9629 
Dorifice  =  0.04819"  ■ 


Dorifce  =  0.07823" 


KL  =  0.9199(V1/V2)2  +  0.1209(V1/V2)  +  0.6892 
R2  =  0.9873 
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Figure  11  -  Correlation  of  HL  Coefficient  for  90  Degree  Sharp-Edge  Orifice 

2 

The  constants  and  R”  values  for  all  L/D’s  are  shown  in  their  respective  Figures. 


Figure  12  -  Correlation  of  HL  Coefficient  for  90  Degree  Sharp-Edge  Orifice 

For  the  L/D  of  2  the  data  where  separation  occurred  tended  to  deviate  from  the  single 
curve  at  higher  velocity  ratios.  It  is  believed  that  this  is  an  artifact  of  the  unstable  nature 
of  the  flow  at  this  condition  and  therefore  it  was  not  considered  in  the  final  correlation. 
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Figure  13  -  Correlation  of  HL  Coefficient  for  90  Degree  Sharp-Edge  Orifice 

The  results  shown  in  Figures  11-13  are  only  for  attached  flow  and  therefore  no  data  are 
shown  for  L/D  =  1 .  Separation  did  occur  at  several  test  conditions  for  the  L/D  of  2  and  5 
orifices  and  this  variable  is  discussed  in  the  Separation  section  of  this  paper. 

INCEPTION  OF  CAVITATION 

As  back  pressure  is  lowered,  the  acceleration  forces  at  the  orifice  entrance  increases, 
thereby  lowering  the  static  pressure.  At  the  point  where  the  onset  of  cavitation  occurs  (~ 
Kcav  =  1.8),  vapor  bubbles  are  observed  fonning  in  the  boundary  layer8  and  15-19  and  Cd 
starts  to  decrease  as  it  transitions  to  full  cavitation  (i.e.  constant  Cc  and  choked  flow), 
Figure  14.  Note  that  the  orifice  cavity  remains  at  the  vapor  pressure  until  the  acceleration 
forces  diminish  sufficiently  for  the  reattachment  process  to  dominate.  After  reattachment 
the  flow  is  similar  to  that  of  the  non-cavitation  characteristics.  As  the  back  pressure  is 
further  lowered  the  acceleration  forces  increase,  and  depending  on  the  orifice  length, 
either  full  cavitation  or  supercavitation  occurs. 

o 

Koivula  noted  that  at  the  inception  of  cavitation  the  fonnation  and  collapse  of  the 
bubbles  resulted  in  both  a  pressure  oscillation  frequency  and  overpressure.  Fie  found 
frequencies  of  8  KHz  at  inception  of  cavitation  but  low  overpressure  fluctuations.  The 
low  overpressure  is  thought  to  be  due  to  the  limited  number  and  size  of  vapor  bubbles  in 
this  regime. 
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Figure  14  -  Illustration  showing  Pressure  Characteristics  for  Inception  of  Cavitation 

The  inception  of  cavitation  is  defined  as  the  point  where  non-cavitation  flow  transitions 
into  cavitation.  This  can  occur  at  a  point  where  Cc  is  immediately  constant  or  transitions 
to  the  point  where  Cc  becomes  constant.  The  length  of  transition  depends  on  the  rapidity 
of  generation  of  bubbles  until  a  vapor  cavity  forms.  Ganippa  noted  that  for  an  orifice 
design  similar  to  that  of  Configuration  #1  the  generation  of  vapor  bubbles  at  inception  of 
cavitation  occurred  completely  around  the  orifice  periphery  (symmetric),  while  for 
Configuration  #2  they  occurred  on  the  left  side  (see  Figure  14)  of  the  entry  at  the  corner 
of  the  orifice  entrance  and  were  observed  only  on  that  side  of  the  orifice  (non- 
symmetric).  In  both  cases  it  appeared  that  the  formation  of  a  cavity  occurred  rapidly  as 
was  observed  in  this  study. 

The  data  show  that  for  some  conditions  there  is  no  bubble  formation  regime  and  the  flow 
goes  directly  from  non-cavitation  to  full  cavitation.  At  other  conditions  there  is  a  clear 
distinction  between  non-cavitation  and  the  Kcav  value  for  full  cavitation  (i.e.  Cc  = 
constant).  This  is  best  shown  by  plotting  the  flowrate  vs  Kcav  and  noting  the  slope 
change  between  the  two  regimes  until  the  flowrate  becomes  constant.  This  change  in 
slope  indicates  the  inception  of  cavitation  and  the  constant  flowrate  indicates  where  full 
cavitation  occurs.  At  low  manifold  pressures  (i.e.  100  psig)  the  orifice  flow  characteristics 
in  this  transition  zone  tend  to  be  erratic  suggesting  instability  of  the  flow.  This  trend 
however  quickly  stabilizes  as  the  manifold  pressure  is  increased.  When  the  flow  is 
erratic  it  is  difficult  to  define  the  inception  point  however  it  is  clear  that  it  is  not  constant 
and  varies  with  velocity  ratio.  At  higher  manifold  pressure  in  some  cases  there  is  a 
transition  zone  and  when  this  occurs  the  inception  point  is  definable  and  is  not  a  function 
of  velocity  ratio.  What  is  causing  this  erratic  behavior  is  not  clear  and  would  require 
visual  observations  of  the  flow  as  well  as  additional  measurements  to  better  understand. 

Understanding  this  phenomenon  is  extremely  important  to  rocket  engine  design  since 
these  erratic  flow  conditions  will  always  occur  during  startup  and  if  the  engine  is  deep 
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throttled  to  a  chamber  pressure  less  than  —  110  psia.  However,  for  this  paper  since  the 
point  of  inception  is  not  easily  determined  when  the  flow  is  erratic  the  inception  point  is 
defined  at  this  time  only  for  well  ordered  flow  conditions. 

The  inception  of  cavitation  for  all  test  conditions  indicated  a  cavitation  inception  at  about 
Kc  av  of  1.8  (Kcav0  5  =  1  .34).  Typical  results  are  shown  in  Figure  15  and  for  this  case  the 
flowrate  was  constant  with  further  lowering  of  the  cavitation  number.  In  some  instances, 
however,  while  the  initiation  of  cavitation  occurred  at  the  same  value  of  Kcav  it  was  not 
abrupt  but  transitioned  to  a  constant  Cc  value.  The  flowrate  also  did  not  become  constant 
until  the  Cc  became  constant.  The  transition  occurred  between  Kcav  of  1 .75  - 1 .8. 


Kcav05 


Figure  15  -  Illustration  of  Inception  of  Cavitation  for  Configuration  #1  and  #2 

The  cavitation  number  inception  for  Configuration  #3  did  not  measurably  change  with 
cross  flow  from  that  found  for  the  other  configurations  although  the  transition  is  less 
pronounced.  This  is  shown  in  Figure  16  where  the  transition  was  between  a  cavitation 
number  of  1.8  to  1.75  (Kcav0'5  between  1.32-1.35). 
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Figure  16  -  Cavitation  Inception  and  Transition  Region  to  Full  Cavitation 


It  is  interesting  that  at  the  point  of  inception  of  cavitation  the  flowrate  is  still  not  constant, 
as  shown  in  Figure  17,  until  the  cavitation  number  reaches  the  point  where  Cc  is  constant 
(Kcav0  5  ~  1.28).  This  finding  is  extremely  important  in  that  it  suggests  that  the  orifice  is 
not  choked  as  previously  hypnotized  by  others  until  full  cavitation  is  reached. 


Figure  17  -  Cavitation  Impacts  on  Flowrate 
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It  should  be  noted  that  Mishra  and  Peles17  and  Yan20  determined  that  the  amount  of 
dissolved  oxygen  in  the  fluid  changed  the  value  of  Kcav  where  the  inception  of  cavitation 
occurred.  The  Mishra  and  Peles  study  for  micro  channels  showed  that  inception  ranged 
from  1.2  to  1.6  as  the  concentration  varied  from  2  to  15  PPM.  They  also  noted  that  the 
transition  in  the  flow  regimes  from  incipient  cavitation  and  then  to  choked  flow  was 
rapid.  The  amount  of  dissolved  gases  in  the  water  used  in  the  tests  of  our  study  was  not 
measured,  however  the  water  used  in  the  tests  was  stored  for  long  periods  of  time  in  air  at 
atmospheric  temperature  and  pressure,  so  it  is  reasonable  to  assume  the  amount  of 
dissolved  gases  was  near  the  equilibrium  value,  which  for  nitrogen  at  298K  is  23  PPM. 
The  equilibrium  value  increases  substantially  as  the  water  is  pressurized.  Care  was  taken 
that  the  water  did  not  spend  long  periods  of  time  in  a  pressurized  state  before  a  test.  In 
addition,  the  piping  configuration  ensured  a  slug  flow  condition,  so  any  increase  in 
dissolved  gases  during  pressurization  would  have  been  limited  to  a  region  near  the  gas  / 
liquid  interface  in  the  tank,  due  to  the  slow  rate  of  diffusion.  Therefore,  it  is  reasonable  to 
assume  that  the  amount  of  dissolved  gases  in  the  water  actually  flowed  through  the 
orifice  was  that  of  the  equilibrium  unpressurized  state.  The  observed  Kcav  value  of  1.8  at 
cavitation  inception  is  therefore  not  inconsistent  with  the  results  of  Mishra  and  Peles. 


For  the  configurations  studied,  within  experimental  accuracy,  the  inception  of  cavitation 

was  found  to  occur  at  the  same  value.  This  was  compared  with  that  predicted  by  the 
2 1 

relationship”  : 


Kincept  =  \.9{\-rlD)1-1^-  (18) 

Rch 

Where  for  our  study  r/D  =  0.003  and  Reh  is  the  “head”  Reynolds  Number  defined  as: 

Re  h  =  Dpi  n yj2g(Px  -  P2)  / p  (19) 

Applying  values  from  this  study  the  Kincept  is  1.88  which  is  close  to  that  experimentally 
determined  and  is  certainly  within  the  accuracy  of  the  data  (Kcav  =  1 .8). 


FULL  CAVITATION  (CONSTANT  FLOWRATE  &  Cc) 

At  cavitation  the  pressure  in  the  vapor  cavity  is  below  the  vapor  pressure  then  rises  as  the 
flow  expands  to  a  pressure  greater  than  the  back  pressure  at  reattachment.  The  flow  then 
slowly  drops  to  the  back  pressure  at  the  orifice  exit.  The  pressure  variation  throughout 
the  orifice  length  for  the  core  and  near  wall  flow  is  illustrated  in  Figure  18.  Full 
cavitation  is  defined  as  the  condition  where  the  flowrate  and  Cc  become  constant  and  the 
attachment  point  has  not  reached  the  orifice  exit. 


18 


Wc,  Pc,  Vc 


Orifice  Flow  Characteristics 


W1,P1,V1  I 


dc 


Wout,  Pout,  Vout 


b 


■  Pv  -  Vapor  Cavity 


Variable 

Reattachment 

Point 


W2,  P2,  V2 
Kcav  <1.3 

Orifice  Static  Pressure  Profile 


Figure  18  -  Static  Pressure  Characteristics  at  Cavitation 

As  the  back  pressure  is  lowered  (P2  approaching  Pv)  the  contraction  coefficient  remains 
constant  but  the  attachment  point  moves  toward  the  orifice  exit4, 15  and  1S.  In  the  limit,  it 
has  been  observed  that  the  attachment  is  at  the  orifice  exit.  During  this  process  the  flow 
at  the  vena-contracta  is  “choked”  and  the  back  pressure  no  longer  influences  the  flow 
parameters  upstream  of  the  vena-contracta.  Since  Cc  is  constant  at  full  cavitation  the  area 
of  the  vena-contracta  also  remains  constant. 

As  the  expansion  process  progresses  droplets  form  in  the  vapor  pocket,  are  re-entrained 
in  the  bulk  flow,  and  new  droplets  formed.  The  amount  of  droplets  in  the  vapor  pocket 
increase  as  the  reattachment  point  approaches  the  orifice  exit. 

If  the  orifice  length  is  insufficient  to  meet  the  length  required  for  reattachment  between 
the  onset  of  cavitation  and  a  cavitation  number  of  1,  then  one  of  two  conditions  occur.  If 
the  flow  is  unable  to  reattach  then  the  flow  will  experience  separation  (hydraulic  flip).  If 
the  flow  does  not  separate,  then  the  unsteady  nature  of  the  flow  at  separation  results  in 
vapor  and  water  exiting  the  orifice.  As  discussed  in  the  Supercavitation  section  of  the 
paper,  a  vapor  cloud  surrounding  the  core  liquid  has  been  observed  and  photographed  by 
other  researchers. 

In  Figure  19  the  results  of  Sato  and  Saito19  for  L/D  of  4.5  illustrates  the  relationship 
between  cavitation  number  and  cavity  extent  in  the  orifice.  It  should  be  noted  that  for 
consistency  Saito ’s  definition  of  cavitation  number  was  converted  to  the  definition  used 


19 


in  this  paper.  These  results  are  extremely  interesting  in  that  they  show  that  the  cavity 
reaches  the  exit  of  the  orifice  at  a  cavitation  number  of  ~1.7.  Sato’s  data  for 
Configuration  #land  L/D  of  2.2  shows  that  full  flow  occurs  at  a  cavitation  number  of  1 .7, 
and  that  decreasing  the  orifice  size  significantly  increases  the  cavitation  number  for  full 
flow.  Also  of  interest  is  that  once  the  cavity  reaches  about  14  of  the  orifice  length  the 
cavity  reattachment  instantly  “jumps”  to  the  orifice  exit.  This  type  of  dramatic  change  is 
mirrored  in  Figure  7  where  Cd  dramatically  decreases  then  becomes  near  constant. 


Figure  19  -  Relationship  between  Kcav  &  Cavity  Extent19 

Full  cavitation  occurs  for  Configurations  #1  and  #2  at  a  Kcav  of  about  1.7-1.75.  The 
results,  as  indicated  in  Figure  20,  clearly  show  that  the  vena  contracta  between 
Configuration  #1  and  #2  are  identical  with  Cc  (0.615).  These  values  are  reasonably  close 
to  Sato’s  results  considering  the  difference  in  the  orifice  diameter. 


Figure  20  -  Comparison  of  Cavitation  Characteristics,  Configurations  #1  and  #2 
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This  suggests  that  while  there  is  a  significant  difference  between  the  cavity 
characteristics,  it  does  not  impact  Cc  (i.e.  Ac).  Note  that  in  all  cases  the  flow  is  attached 
and  cavitation  remains  to  at  least  a  value  of  Kcav  =  0.996.  It  should  be  noted  that  the 
minimum  back  pressure  was  0.042  psig.  For  this  orifice  design  the  expected  Cc  is 
determined  fronT: 


C  = 


0.61 12 


—  11.4 


D 


(20) 


For  the  orifices  used  in  this  study  r/D  is  <  .003  resulting  in  a  value  of  Cc  of  0.615. 
Configuration  #1  experimentally  determined  value  compares  with  that  determined  from 
Equation  20.  Also,  the  experimentally  determined  Cd  of  0.76  compares  well  with  the 
data  of  Licharowicz  (0.75). 

As  indicated  in  Figure  21  the  relationship  between  Cd  and  Cc,  Kcav  is  still  linear  even  with 
cross  flow. 


Figure  21  -  Constant  Contraction  Coefficient  in  a  Sharp  Edge  90  degree  Orifice 

Results  shown  in  Figure  22  validate  that  when  Cc  is  constant  the  flowrate  is  also  constant. 
This  is  similar  to  that  of  a  cavitating  venturi. 
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Figure  22  -  Illustration  of  Constant  Flowrate  in  the  Cavitation  Regime 

Inspection  of  Equation  13  shows  that  the  data  analysis  for  Cc  is  at  constant  A1/A2,  V 1/V2, 
and  L/D.  Therefore,  the  calculated  Cc  may  be  independently  a  function  of  these  variables: 


Cc=f(A1/A2,V1/V2,L/D)  (21) 

From  a  physical  basis  the  area  ratio  should  only  impact  the  flow  contraction  and  L/D 
would  only  impact  the  expansion,  turbulence,  and  friction  losses  flow  downstream  of  the 
vena-contracta.  The  independent  impact  of  velocity  ratio  should  also  impact  head  loss. 
These  losses  all  impact  the  overall  head  loss  and  if  significant  would  result  in  an 
independent  impact  on  Cc. 

To  test  these  variables  the  data  are  plotted  in  Figures  23,  24,  and  25  as  function  of 
velocity  ratio  and  area  ratio  for  a  specific  L/D.  The  results  show  that  at  L/D  of  10,  Cc  is 
not  a  function  of  area  ratio  although  as  the  L/D  decreases  there  is  an  indication  of  a  minor 
impact.  This  is  the  result  of  the  flow  not  following  the  constant  Cc  characteristic  at  the 
lower  value  of  L/D  suggesting  that  the  flow  is  being  impacted  by  other  processes  not 
included  in  the  original  derivation.  If  the  L/D  characteristics  are  compared  there  is  also 
indication  of  an  independent  impact  although  the  variation  is  within  the  overall 
experimental  error.  Note,  also  that  Cc  tends  to  decrease  as  the  velocity  ratio  is  increased 
beyond  ~  0.45.  For  velocity  ratios  greater  than  about  0.5  represents  the  data  obtained  at 
an  upstream  pressure  of  100  psi,  which  again  is  in  the  flow  conditions  where  the  flow 
was  not  well  behaved.  It  is  possible  that  these  variables  have  an  independent  impact, 
however  sufficient  measurements  were  not  taken  to  identify  what  is  causing  the  flow 
behavior  in  the  low  upstream  pressure  and  short  L/D. 
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Figure  23  -  Impact  of  Area  Ratio  and  Velocity  Ratio  on  Cc  for  L/D  =  10 
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Figure  24  -  Impact  of  Area  Ratio  and  Velocity  Ratio  on  Cc  for  L/D  =  5 
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Figure  25  -  Impact  of  Area  Ratio  and  Velocity  Ratio  on  Cc  for  L/D  =  2 

For  well  behaved  flow  conditions  Equation  13  provides  excellent  prediction  of  the 
contraction  coefficient.  Flowever,  for  conditions  where  the  flow  may  be  less  well 
behaved  Figures  23-25  provide  an  estimate  of  the  additional  impacts.  The  development 
of  a  complete  correlation  requires  additional  data  with  emphasis  in  the  low  upstream 
pressure  range  and  L/D  less  than  5. 

SUPERCAVITATION 

Chaves  et  al15  defines  supercavitation  to  describe  when  the  cavities  extended  past  the  exit 
of  the  nozzle  but  the  orifice  does  not  experience  separation  and  the  spreading  jet  is 
surrounded  by  a  vapor  cloud  as  indicated  in  Figure  26.  Chaves  also  noted  that 
supercavitation  was  accompanied  by  a  dramatic  increase  in  the  spray  angle  produced  by 
the  nozzle.  Visual  observations7  of  cavitation  flows  in  nozzles  have,  in  addition,  shown  a 
transient  behavior  caused  by  quasi-periodical  re-entrant  jet  motion  in  the  vena-contracta 
region,  and  random  collapse  of  the  cavitating  bubbles  at  the  nozzle  exit. 
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Figure  26  -  Illustration  of  Supercavitation  Jet  Characteristics 

The  two-phase  bubble  nature  of  cavitation  flows  has  been  confirmed  by  measurement  and 
visualization7'  8’  15  and  19.  At  this  condition  the  flow  is  no  longer  choked  and  Cd 
experiences  either  a  dramatic  step  decrease  followed  by  more  gradual  lowering  (as  the 
back  pressure  is  further  lowered)  until  the  minimum  is  reached  or  complete  separation 
allowing  the  back  pressure  to  surround  the  jet  from  the  orifice  exit  to  the  vena-contracta. 

o 

Koivula’s  results  in  this  regime  show  that,  due  to  the  implosion  of  vapor  in  the  cavity 
which  produces  shock  waves,  the  frequency  of  oscillations  decrease  from  8  KHz  at 
inception  of  cavitation  to  -700  Hz  at  supercavitation  and  the  overpressure  increases  to  as 
high  as  H^90  psi. 

In  tests  utilizing  Configuration  #2  with  the  L/D  of  5  orifice,  supercavitation  occurred  at 
Kcav  of  about  1 .6  (Kcav°5  =  1  .25;  see  Figure  7)  and  at  this  point  the  discharge  coefficient 
dropped  dramatically.  Upon  further  reduction  the  decrease  in  Cd  tended  to  reach  a 
constant  value  equal  to  that  of  separation.  Ahn  et  al5  also  noted  that  for  an  L/D  of  5  and 
orifice  size  of  0.04”  as  the  cavitation  number  is  increased  the  flow  transitions  very 
quickly  from  separated  flow  to  cavitating  flow  similar  to  that  shown  in  Figure  7. 

The  data  in  Figure  27  compares  the  Configuration  #1  and  #2  cavitation  characteristics. 
For  Configuration  #2  Cc  is  constant  in  the  cavitation  regime  while  for  Configuration  #1 
supercavitation  occurs  as  noted  by  the  variation  from  the  linear  curve. 


25 


Figure  27  -  Illustration  Showing  the  Gradual  Onset  of  Supercavitation 

The  experimental  data  confirm  that  the  flow  is  not  constant  in  the  supercavitation  regime 
suggesting  that  the  cavity  is  no  longer  choked  as  shown  in  Figure  28.  This  is  similar  to 
the  onset  of  cavitation  where  the  flowrate  varies  but  returns  to  the  prior  value 
corresponding  to  that  at  cavitation. 


Figure  28  -  Variable  Flowrate  in  the  Supercavitation  Regime 

One  explanation  for  the  variation  in  flowrate  is  that  the  pressure  oscillations  at  the  orifice 
exit  caused  by  the  unsteady  re-entrainment  process  results  in  the  cavitation  number 
varying  sufficiently  that  the  flow  varies  between  cavitation  inception,  cavitation,  or  the 
non-cavitation  regimes.  For  our  experiments  no  high  speed  instrumentation  was  utilized 
in  these  experiments  to  determine  the  magnitude  and/or  frequency  of  the  oscillations. 

o 

However,  considering  the  magnitude  of  overpressure  measured  by  Koivula  ,  the  flowrate 
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could  easily  vary  from  0  to  choked  flow  to  separation  at  very  high  frequencies.  The 
mean  flowrate  of  course  would  be  lower  than  that  of  choked  flow  similar  to  that  shown  in 
Figure  18.  Also  the  high/low  overpressure  combined  with  the  very  high  frequencies 
could  account  for  the  vapor  leaving  the  orifice  and  not  allowing  the  back  pressure  to  enter 
the  orifice.  This  is  an  area  that  requires  additional  study  and  experimentation. 

The  cross  flow  results  tended  to  be  less  impacted  by  supercavitation  then  the  other 
configurations.  In  Figure  29  the  Configuration  #3  data  where  velocity  ratio  was  varied 
show  that  supercavitation  did  not  occur.  This  is  contrasted  with  Figure  27,  Configuration 
#2. 


Figure  29  -  Configuration  #3  Cross  Flow  Impacts  on  Supercavitation 
SEPARATION  (HYDRAULIC  FLIP) 

Separation  occurs  when  the  acceleration  forces  at  the  entrance  to  a  sharp-edge  orifice  are 
sufficient  to  cause  the  flow  to  fully  detach  from  the  wall  and  the  orifice  L/D  is 
insufficient  to  allow  for  reattachment.  An  illustration  of  separation  is  provided  in  Figure 
30. 

For  the  case  where  the  flow  fully  separates  the  Cd  will  immediately  decrease  to  the 
minimum  value  and  remain  there  as  the  back  pressure  is  lowered.  An  example  of  the 
dramatic  drop  in  Cd  at  separation  then  a  near  constant  Cd  as  the  cavitation  number 
approaches  1  is  shown  in  Figure  31.  As  stated  above  the  initiation  occurs  at  Kcav(l'5  -  1.25 
(cavitation  umber  of  ~1.6).  The  variation  in  Cc  in  the  cavitation  regime  is  due  to  the 
velocity  ratio  impacts  as  discussed  above. 
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Figure  30  -  Illustration  of  Separation  (Hydraulic  Flip) 


Figure  31  -  Flow  Characteristics  Illustrating  Separation 

The  variation  in  orifice  flowrate  for  the  60  ft/sec  data  is  shown  in  Figure  32.  These 
results  show  that  the  flow  is  constant  in  the  cavitation  regime  then  experiences  a  step 
change  in  flowrate  at  the  onset  of  supercavitation,  and  then  continues  to  vary  as  the 
cavitation  number  approaches  1 .  Further  note  the  large  shift  in  flowrate  rather  than  the 
gradual  change  that  is  observed  in  the  supercavitation  regime. 
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Figure  32  -  Illustration  of  Flowrate  Variation  at  Separation 

2 1 

A  relationship  for  the  critical  cavitation  number  for  hydraulic  flip  is: 


Kcrit  =  1  + 


_ 1 _ 

(1 +  LI  4D)(1  +  2000  /  Re  h)e70r,D 


(22) 


The  results  from  this  study  are  compared  with  the  predictions  from  Equation  22  in  Figure 
33.  Note  the  excellent  agreement.  The  L/D  of  10  for  the  0.07823  inch  orifice  does  not  fit 
the  correlation  since  it  did  not  experience  separation  within  the  range  of  this  study. 
However,  the  smaller  orifice  did  experience  separation  and  correlates  well  with  the 
prediction.  Equation  22  is  stated  to  be  valid  only  for  short  L/D  orifices  ,  and  these 
results  suggest  that  the  limit  is  about  L/D  =  10. 


Figure  33  -  Comparison  of  measured  Separation  with  Equation  22  Correlation21 
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The  contraction  coefficient  at  separation  is  shown  in  Figure  34.  The  correlation  should 
be  used  with  caution  since  the  range  of  the  parameters  is  limited.  The  resulting  equation 
is: 

C  =-0.1054^  +  0.6158  (23) 

V2 


Figure  34  -  Contraction  Coefficient  vs  Velocity  Ratio  at  Separation 

CONCLUSIONS 

The  results  presented  clearly  show  the  complexity  of  possible  flow  conditions  for  sharp- 
edge  orifices  in  cross  flow.  For  example,  an  orifice  can  initially  flow  separated  then 
either  (1)  quickly  attach  in  the  cavitation  regime,  (2)  attach  in  the  non-cavitation  regime, 
or  (3)  remain  separated  regardless  of  flow  conditions.  The  ability  to  select  which 
scenario  the  orifice  will  be  operating  in  is  the  challenge.  Once  this  is  selected  then  all 
pertinent  flow  parameters  can  be  defined. 

The  most  significant  contribution  of  this  study  is  the  extension  of  the  linear  cavitation 
model  verified  in  numerous  studies  to  include  both  cavitation  and  non-cavitation  flows 
with  3  manifold  fluid  feed  configurations.  The  supporting  relationship  provides  design 
engineers  with  the  ability  to  more  knowledgably  predict  flow  regimes  as  well  and  flow 
variables  for  determination  of  CCj  Cd  and  KL. . 

In  the  non-cavitation  regime  the  results  of  this  study  show  that  the  head  loss  can  be  as 
high  as  8  times  that  for  no  cross  flow  depending  on  the  magnitude  of  the  cross  flow.  This 
emphasizes  to  the  design  engineer  the  need  to  keep  the  cross  flow  velocity  low  to 
minimize  pressure  requirements.  In  addition  the  data  indicate  that  the  major  head  loss  in 
the  non-cavitation  regime  occurs  upstream  of  the  vena-contracta. 

The  inception  of  cavitation  and  the  transition  to  full  cavitation  has  been  identified. 
However,  additional  studies  are  required  to  quantify  the  variables  impacting  the  transition 
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from  initiation  of  bubbles  to  the  full  cavity  formation.  Visual  studies  have  also  identified 
the  transition  and  have  shown  that  it  can  transition  to  full  cavity  very  rapidly.  No 
correlations  have  been  offered  to  date. 

In  the  cavitation  regime  the  validity  of  the  first  order  relationship  was  shown  to  be  related 
to  the  fact  that  the  head  loss  from  the  manifold  to  the  vena-contracta  is  small.  Cross  flow 
impacts  Cc  revealing  the  significance  of  cross  flow  velocity  on  the  vena-contracta 
formation  process. 

This  study  has  determined  that  the  flow  is  not  choked  in  the  supercavitation  regime.  It  is 
postulated  that  this  is  related  to  the  unsteady  nature  of  the  flow  (frequency  and 
overpressure).  In  addition  the  onset  of  supercavitation  is  not  defined  sufficiently  so  that 
with  certainty  a  prediction  that  it  will  occur  and  further,  when  it  occurs,  will  it 
immediately  cause  separation  or  quasi-attachment  until  separation  occurs.  The  results  of 
this  study  can  only  determine  that  once  supercavitation  occurs  it  CAN  cause  immediate 
separation  or  quasi  separation  as  back  pressure  is  further  lowered.  Additional  studies  are 
required  to  better  define  where  supercavitation  will  occur  by  defining  the  conditions  that 
cause  the  cavity  to  grow  and  suddenly  jump  to  the  orifice  exit  and  then  the  change  in  the 
processes  occurring  at  the  exit  leading  to  separation.  Of  particular  interest  to  rocket 
engine  injector  design  is  the  potential  for  the  unsteady  nature  of  supercavitation  to  trigger 
combustion  instability  in  the  combustion  chamber  especially  during  startup. 

The  impact  of  L/D  for  sharp-edge  90  degree  orifices  has  been  defined  and  the  data  was 
correlated  with  an  existing  correlation  predicting  separation  with  short  L/D  orifices.  This 
correlation,  however,  can  only  predict  that  separation  is  likely  to  occur. 

NOMENCLATURE 


SYMBOLS 


2 

A  Area  in";  Constant  in  Equation 

B  Constant  in  Equation 

C  Constant  in  Equation 

Cc  Contraction  Coefficient,  A2/A1 

Cd  Discharge  Coefficient 

D  Orifice  Diameter,  in 

Hl  Head  Loss,  lb/in“ 

K  Loess  Coefficient 

Kcav  Cavitation  Number 

L  Orifice  Length,  in 

P  Pressure,  psi,  psia 

Re  Reynolds  Number 

r  Orifice  Entry  Radius,  in 

V  Velocity,  ft/sec 

W  Flow  Rate,  lb/sec 
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GREEK 


p  Liquid  Density,  lb/ft 

SUBSCRIPTS 

1  Manifold 

2  Orifice  Exit 

a  Contraction 

c  Contraction 

crit  Critical  at  Separation 
h  Head 

incept  Cavitation  Inception 
v  vapor 

out  Exiting  Manifold 

turb  Turbulence 

f  Friction 
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